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Abstract. We investigate Dilcher and Stolarsky's polynomial analogue of the Stern 
diatomic sequence. Basic information is obtained concerning the distribution of their 
zeros in the plane. Also, uncountably many subsequences are found which each converge 
to a unique analytic function on the open unit disk. We thus generalize a result of Dilcher 
and Stolarsky from their second paper on the subject. 



1. Introduction 

The Stern sequence {a{n)}n>Q, also known as the Stern diatomic sequence (or Stern's 
diatomic series), is an integer sequence defined recursively by a(0) = 0, a(l) = 1, and 

a{2n) — a{n), 

a{2n + 1) = a{n) + a{n + 1). 

Dilcher and Stolarsky recently introduced a polynomial analogue of this sequence and ex- 
plored some of its properties [SJ [B] (see [S] for background information and references for 
the classical Stern sequence). This new sequence {a(n;z)}„>o is defined by a(0;z) = 0, 
a(l; z) = 1, and 

a(2n: z) — a(n: z^), 
(1) ^ ' ; V , 

a{2n + 1; 2) = z a{n; z'^) + a{n + 1; z^). 

These polynomials have only and 1 as coefficients, and a{n: 1) = a(n). See Table [2] for a 
list of the first few polynomials in the sequence. 

As nothing has yet been said about the zeros of these polynomials, we discuss the basics 
of their distribution in the plane in Section [2j 

The heart of this paper lies in Section [3j where we switch to a more structural point of 
view to show that subsequences of ([I]) of a certain form converge to analytic functions on the 
open unit disk. In fact, there are uncountably many such limit functions. This generalizes 
a result of Dilcher and Stolarsky in [B] where two particular convergent subsequences of ([T]) 
were found. The transcendence of the limit functions of these subsequences has been of 
interest: Coons 4 showed that the limit functions are transcendental over Q{z), the field 
of rational functions over Q, and Adamczewski [T] further showed that the functions take 
transcendental values at every algebraic number in the punctured unit disk. Also, Fibonacci- 
like quadratic relations were derived in |n| for the elements of these subsequences — an effort 
continued by Vsemirnov in |12| . We will find these particular subsequences and functions 
among those introduced in this paper in Section |4] 

We take a moment here to note that a different polynomial analogue to the Stern sequence 
was given by Klavzar et. al. in [9]. A similar sequence to this was defined by Alkauskas [2] 
in relation to p-continued fractions. 
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2. Zeros 

In this section we show that the zeros of the Stern polynomials cluster uniformly near 
the unit circle in the sense of Weyl. That is, we can prove the following result. 

Proposition 2.1. Let i'n{p) he the number of zeros of a{n; z) in the annulus 
I — p < l^l < 1/(1 — p), and let Vn{Q, <P) he the number of zeros of a{n; z) in the sector 
9 < arg z < (p. For fixed p, 6, and (j) satisfying < p < 1 and < 9 < (f> < 27r, 



/ 31ogn-log8 ^ iyn{0,(l)) _ - g ^ / 31ogn-log8 
V n deg a{n;z) 2tt \ n 

and 

< 1 _ '^"(P) < 2 3 log n - log 8 
~ dega(n;z) p n 
for all n > which are not powers of 2. 

To prove this we require three results, the first due to Dilcher and Stolarsky [5], the 
second due to Erdos and Turan |7j, and the third due to Hughes and Nikeghbali 0. 

Proposition 2.2 (Dilcher and Stolarsky). Let e{n) = ord2(n), the highest power of 2 
dividing n. Then for n > 1, 

, , n-2^(") 
dega(n;z) = . 



Proposition 2.3 (Erdos and Turan). For a polynomial 

N 

P{z)=Y,a,z= 

with aoa^ 7^ 0, define 

^ 1 1 

L{P) = log^lajl - -log|ao| - -log|aAr|. 

If vp{9, <j)) is the number of zeros of P{z) in the sector 9 < arg z < <f>, then, for Q < 9 < (f> < 
2-K, 

Proposition 2.4 (Hughes and Nikeghbali). For a polynomial 

N 

P(z)=^«,z^' 



with a^aN 7^ 0, let L{P) be as defined in Proposition 2.3 If i>p{p) is the number of zeros 



of P{z) in the annulus 1 — p < |z| < 1/(1 — p), then, for < p < 1, 

up{p)\ ^ 2 L{P) 



N - p N 



We may now prove Proposition |2.1[ 
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Proof of Proposition 2.1 In light of Propositions 2.3 and 2.4 the resuhs will follow from 
estimating the quantity L{a{n; z))/ deg a(n; z). 

Let e{n) be the largest power of 2 dividing n, and let d{n) — n/2'^(") be the odd component 
of n. In these terms, d{n) = 1 is necessary and sufficient for a(n;z) — 1, so we will only 
consider n such that d{n) > 3. With the aid of Proposition |2.2| we calculate 

L{a{n; z)) < log deg a{n; z) 
< log n — log 2 



and 

to get 



deg a(n;z) > - 



L{a{n;z)) 3 log n — log 8 
deg a(n; z) n 



This concludes the proof. □ 
We conclude this section with a small remark about real zeros. 

It follows from Proposition |2.2| that the only Stern polynomials of odd degree are those of 
the form a(4n + 3; z). Now, a(4rt + 3; — 1) — —a{n) < and a(4n + 3; 0) = 1, so a(4n + 3; x) 
(where x is real) has a zero in the interval (—1,0). Further, using ([T]) it is not difficult to 
show that on (— oo, —1) we have a{4n + 3; x) < 0, on (0, oo) we have a(4n + 3; x) > 0, and 
on (—1,0) we calculate 

-^a(4n + 3;x) = (3x'^ + ix^) a'(n; x^) + (1 + ix^ + 4a;^) a'(n + 1; x'^) > 0. 
ax 

Hence this real zero of a(4n + 3; z) is unique. 

3. Convergent subsequences 

Given the definition of this polynomial sequence it is natural to look at subsequences of 
the form {a(2"+in + ^; z)}„>o for fixed integers m > and < £ < 2™+i. To facilitate this 
we will enlist the help of binary sequences (that is, sequences whose elements are only O's 
and I's). We will denote a finite sequence by {bk}o<k<m and an infinite sequence by {bk}- 

Suppose we have a binary sequence {6fe}o<fc<m- With this sequence we may associate 
the polynomial sequence 

I / 

n>0 

For fixed to, every polynomial sequence of the form {a(2™+^ri + i; z)}n>Q is associated with 
exactly one finite binary sequence. 

But we are more interested in what occurs when we instead leave n fixed, begin with an 
infinite binary sequence {&fe}, and allow to to vary, producing the sequences 




m>0 

These sequences are unique in that as to increases the heads of the polynomials stabilize, 
so that they gradually build a formal Maclaurin series. The series built turns out to be 
independent of n, thus we identify it with the binary sequence {bk}- 
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Proposition 3.1. Given a binary sequence {bk}, the sequences of polynomials 
a| 2'"+^i 2'bj ■ z\ \ for all fixed n > 

converge to the same analytic function f{bk} '^'^ open unit disk. 

The first few terms of one of the subsequences from Proposition 3.1 are shown in Table [l] 
In the context of Hurwitz's Theorem (see e.g. pQQ p. 4]), the resuh is illustrated graphically 
in Figure [l] which displays the convergence of the roots of the polynomials. 
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+ 1 + 
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z6 
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+ + 


zU 


- Zl" + Zl2 
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+ + 


zU 


- z^ + z^^ 4 


-z20 


+ z24 
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+ z^ + 


zU 


- z^ + z^^ 4 


-z20 


+ z24- 




h z52 + ^56 




5 
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+ z2 + 


zU 


- z** + z^^ 4 


-z20 


+ z24_ 


f z32_ 


h + _|_ 


z40 + _^ ^84 






+z** 4 




+ z98 + zio" + 


^104 








6 


1 


+ z^ + 


zU 


- z** + z^* 4 


-z20 


+ z24- 


f z32_ 


h z34 + + 


z40 + + 






+Z684 


-Z72 


+ z"6 + ^148 4. 


^152^ 


^160, 




+ zl68 + zl92 






+ zl94 


f ^196 + ^200 













Table 1. The first few elements in the sequence from Proposition |3. 1 1 with 
n = 3 and {fefc} = {0, 1, 0, 0, 1, 0, 0}. 



The foundation of the proof of Proposition |3.1| is laid by the following lemma, whose 
proof is a straightforward exercise in mathematical induction. 

Lemma 3.1. Let {6fe} be a binary sequence. Define the two polynomial sequences {sk} and 
{tk} by s_i = 1, = 0, and 

Sfe = (1 - bk)sk-i + z^ [bk{sk-i - tk-i) + tk-i], 
tk = bkSk-i + tk^i- 



The 



o 2^+^n +^ 2^bj ; z = s„ a(n; z^""^'^ + i„ a(n + 1; 



Rearranging the formula in Lemma 3.1 we get 



(2) 

where 
(3) 



a 2^+in +^ 2^6^ ; z A„ 



0=0 



a|n; z^ ^ ^ — a^n + 1; z^ ^ ^ 



The following calculation will be useful. 
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Figure 1. Zeros of elements of the sequence from Proposition 3.1 with 
{bk} = {1, 0, 1, 0, 1, 0, 1, 0, 1}. Each plot fixes m and shows the zeros of the 
elements with n = 5,6, 7, 8, 9. The index m increases from in the top left 
to 8 in the bottom right. 



Lemma 3.2. The degree of the lowest power of z in Am is at least 2™+^ + X^jLo ^"'^i ^^^"i^ 
n> 1. 

Proof. Suppose that 1 occurs at least once in the sequence {fefc}o<fc<m and let kg be the 
index of its first occurrence. We then have Skg-i = 1 and tko-i = 0, so that Sko = z^''"- If 
1 occurs at least twice in {bk}o<k<m, let ki be the index of its second occurrence. Then 

fcl— fco — 1 

z 

3=0 
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which tells us that the degree of the lowest power of z in s^^ is 2^° + 2'^^. Continuing in 
this way we see that, if {fcg, fci, . . . , kp\ is the complete set of indices of occurrences of 1 in 
{^fc}o<fe<m) the degree of the lowest power of z in Sm is 

p m 
j=0 3=0 

The result follows. □ 



This, combined with Proposition 2.2 will allow us to prove Proposition 3.1 



Proof of Proposition \37l\ We will first prove the result for ?i = 1. If 1 never occurs in the 
sequence {bk}, then the "polynomials" in question are just the constant function with value 
1, so suppose it occurs at least once. In the notation of Proposition 2^ we have 

ef2'"+i+^^2^6j = ko, 



where kQ is the index of the first occurrence of 1 in {bk}- Proposition 2.2 then tehs us that 

(rn \ m 

2^+1 2^6,. 2™+ J2 
3=0 J 3=ko + l 

Note that 

^" ^ ^m— 1 ^ [^m(^m— 1 ^m— l) ^m— l]- 

From Q we have that the degree of Sm-i + ^m-i is strictly less than 2™, so that the limit 
of Sm + as TO — ^ cx) is a formal power series ^ ctjZ^ with ao = 1 and aj — or aj = 1 for 
all J > 0. By comparison with the geometric series ^ z^ this series converges on the open 
unit disk to an analytic function f{b^}- 

Now, dropping any restrictions on the sequence {6fc}, we will consider the case where 
n > 1. Equation Q implies that 

deg(s„,+i„0 < 2™+i. 



According to Lemma 3.2 the value on the right-hand side of the above inequality is at most 
the degree of the lowest power of z in the polynomial Am defined in ^ and ([s]). Since the 
constant term of a{n + 1; z^ ) is one we can now write 

m \ 

2^bj ; z\ = s,n +tm + higher order terms, 

3=0 ) 

where the higher order terms are o|z^"'^^^ as z — > 0. Hence, for any n > 0, 
lim af 2"+in + V Vhj ; z | = lim (s„, + t^) = f{t,} 

m— foo I — ^ / m— 5-00 ^ ^ 

on the open unit disk. □ 
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It follows from the above proof and Lemma 3.1 that no two distict binary sequences may 
produce the same power series. 

We remark that the zeros of the polynomials seem to converge outside of the unit circle 
as well, but that the series in question diverge there. Carlitz's study [3] of the reciprocals 
of the Stern polynomials may be of some use in exploring this phenomenon. 



4. The subsequences of Dilcher and Stolarsky 
Define the integer sequence {^„} by 

= ^ (2" -(-!)"), 

and the polynomials {(/>«} by 

<l>n{q) = a(Cn;g)- 

These polynomials were studied by Dilcher and Stolarsky in [6] (we continue their use of q 
as the variable). In particular, they defined the limit functions 

F{q) = lim (j)2n{q), 

n^oo 

G{q) = lim (j>2„+i{q). 

n— ^oo 

Let us find these sequences among those considered in the present paper. 
Let {(3n} be the binary sequence {1, 0, 1,0,.. .}. Then 



2ri-3 

-)2n-2 I 



i=o j=o 



I (4" - 1) 



SO that 

S2n-3 + i2n-3 



where s„ and t„ are defined as in Lemma 3.1 for the binary sequence {/?„}. From this we 
see that F — f^^^y 

Let {jn} be the binary sequence defined by 70 = 1 and 7„ = for n > 1, so that 

{jn} — {1, 1, 0, 1,0,.. .}. We proceed similarly, finding 



2ji-2 n 

2--! + ^ 2^7, = 1+2 E4^' 
j=o j=i 

= - (22"+i + 1) 
3 

SO that 

S2n~2 + t2n-2 = 4'2n+l; 



where s„ and t„ are defined as in Lemma 3.1 for the binary sequence {jn}- Thus we have 

G = /{7„}- 
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The study of these subsequences was motivated by the fact that (/)„(!) — Fn, the n**^ 
Fibonacci number [TUl. Inspired by the recurrence F„ = F„_i +F„_2, Dilcher and Stolarsky 
derived the relations 

4'2niq) = (l>2n^l{q^) + Q (f>2n-2{q^) 

and 

4'2n+l{q) = q(f>2n{q'^) + (t>2n-l{q^), 

and, drawing from the identity Fn-iFn+i — F^^ = (—1)", 

02«+l('?)?!'2n-l(g^) - q(l>2n{q)(l>2n{q^) = 1 

and 

<^2«+l(g)'?!'2n+l(9^) - g02n+2(g)02n(g^) = 1- 

Taking the hmit as n goes to infinity immediately yields interesting functional relations 
between F and G, but perhaps it is more crucial that F and G satisfy 



Fiq)Giq' ") - q^'-Giq)F{q^ ") = ^2niq), 
G{q)G{q''"'") - q^^"-^F{q)F{q^""'') = cj^^n^q), 



for all n > 1 jnHHj. The sequence {4>n} is unique in its relation to the Fibonacci sequence, 
but we suspect that {4>2n} and {(t>2n+i} are are not the only subsequences introduced in 
Proposition 3.1 which satisfy these types of identities. 

For example, consider the nonzero binary sequence with the simplest structure, {6n} = 
{1,1,1,...}. If we define the integer sequence {i^„} by 



n-2 



3=0 

and the polynomial sequence {'ipn} by 

tpniq) = a{vri]q), 

then we have 

(5) = 2Vn-l + l = iyn-l+ 2""^ 

and, from Q, 

(6) ^„(g) =gV^„_i(<z2) + l. 
Induction will show that 

(7) il^niq) + q^"^ = ipn+i{q), 

which, when combined with ([6]), yields 

tpn{q)tpn{q'^) - = 

Unfortunately, no functional relation can be obtained from this by letting n ^ oo. On the 
other hand, we also have from ^ the identity 

(8) Vn(9)-g'^„-2(9*) = i + g-V^2(g), 

and if we let n — >■ oo we deduce that 

(9) /{^„}(9)-9V{^„}(/)-^2(g). 
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Equation Q, combined with a simple induction argument using ([s]), Q, and Q, yields the 
following result. 



Proposition 4.1. 

(10) 

for all 2 < 2m < n. 

Holding m fixed in (10) and letting n — > oo gives a generalization of ([9|. 
Corollary 4.1. 

/{<5„}(9) - = "02™ (g) 

for all m > 1 . 
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a(n 










n 


a{n 














1 


1 










17 


1 + 


z + 


z'- 


Fz4 + 


z« 






2 


1 










18 


1 + 


z2- 


f Z4 


+ z^ 








3 


1 + 


z 








19 


1 + 


z + 


z3- 


l-z4 + 


z^ + 


z* 4 


^z9 


4 


1 










20 


1 + 


z4- 


f z* 










5 


1 + 


z + 








21 


1 + 


z + 


z2- 


l-z^ + 


z6 + 


z* 4 


- z^ + zio 


6 


1 + 










22 


1 + 


Z2- 


f z^ 


+ z^ - 


f zio 






7 

8 


1 + 
1 


z + 








23 
24 


1 + 
1 + 


z 4- 

z8 


z3- 


\-z'^ + 


z^ + 


z9 4 


-z" 


9 


1 + 


z + z^-\- 


-z* 




25 


1 + 


z 4- 


Z2- 


^z^ + 


z9 + 


zio 


f zi2 


10 


1 + 










26 


1 + 


z2- 


f Z4 


+ zio 


+ Zl2 




11 


1 + 


z + 


zU 


-z4^ 


hz5 


27 


1 + 


Z + 


z3- 


Fz4 + 


z^ + 


z" 




12 


1 + 










28 


1 + 


z4- 


f Zl2 








13 


1 + 


z + 


zU 


-Z5^ 


h z^ 


29 


1 + 


Z + 


z2- 


l-z^ + 


z6 + 


^13 


f zi4 


14 


1 + 










30 


1 + 


z2- 


f Z^ 


+ z" 








15 


1 + 


z + 




-z^ 




31 


1 + 


Z + 


z3^ 


l-z'^ + 


^15 






16 


1 










32 


1 















Table 2. The polynomials a(n; z) for ti = 1, 2, . . . , 32. 
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